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Abstract 

We computed the Bk parameter on the lattice by using a non-perturbatively 
improved Wilson action. From our quenched simulation we obtain B^ s (2GeV) = 
0.73(7) lo'oi- To gain some insight in the systematic errors due to the use of the 
quenched approximation, we also made the unquenched computation with Wilson 
fermions. We do not observe any significant deviation with respect to the quenched 
result. 
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1 Preliminaries: Motivation and Status 



The experimental value of the parameter ek, which is the measure of indirect CP violation 
in the kaon system, has been accurately determined [l|]. After combining |£^ xp '' ) | with its 
theoretical expression derived in the Standard Model (SM), one obtains a hyperbola in the 
(p, fj) plane g g, 

fjA 2 B K ( 1.11(5) • A 2 (l - p) + 0.31(5) ) = 0.226 , (1) 

which, together with two equations involving 5-mesons, represents one of the main con- 
straints onto the shape of the CKM unitarity triangle. In eq. ([I]), the quantities A = 0.22, 
A, fj , p are the standard Wolfenstein parameters of the CKM matrix elements 0, |5|, while 
Bk is the so-called bag parameter which encodes the low energy QCD contribution to the 
K° — K° mixing amplitude. The main uncertainty in eq. (|l|) comes precisely from Bk (to 
which we focus in this letter) and to less extent to the parameter A which is known from 
\V cb \=A\ 2 . 

We first remind the reader the definition of Bk- As introduced in ref. ||, the parameter 
Bk measures a mismatch between the hadronic matrix element responsible for the K° — K° 
mixing and its value obtained by applying the vacuum saturation approximation (VSA), 
i.e. 

(K°\Q(p)\K°) = ^f K m 2 K B K (p), (2) 

where p indicates the scale dependence of the operator Q = 37^(1 — S7^(l — 75)0?, 
in some renormalization scheme. The QCD anomalous dimension of this operator, in the 
RI-MOM and MS schemes, is known from perturbation theory to NLO accuracy |7|]- |p.U[| . 
To that precision, Bk{p) can be converted to the renormalization group invariant (RGI) 
form, 

B K = a s (p)~ M (l + B K (P) , (3) 

which is the one that is needed in eq. ([I]). Besides the scheme independent 70 = 4 and /5 = 
11 — 2nf/3, (rif is the number of active flavours), the NLO coefficient J is renormalization 
scheme dependent and in the two schemes which we use in this letter, it is given by J7j-[10 



17397 - 2070n/ + 104n^ 
-Tri-mom — 8 log 2 — — 



6(33 - 2n f ) 



2 



13095 - 1626 ti/ + 8 n) 
Jm = 6(33-2 % ) 2 • (4) 

The main task is to compute Bk{p) in some scheme, by using a non-perturbative technique, 
and then to convert the result to the RGI form, B K - We choose lattice QCD and the RI- 
MOM scheme. 
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Before discussing our computation, let us briefly recapitulate the actual situation con- 
cerning the value of B^, which at the same time gives a good motivation to compute this 
quantity. 

Historically, the three distinct approximations applied to this problem lead to different 
results: VSA gave B K = 1 @, PCAC and SU(3) symmetry implied B K = 1/3 Q, and 
the limit of the large number of colours (N c — > oo) yielded B K = 3/4 |12|]. [] Such a 
large spread of values prompted many physicists to develop various models which typically 
combine the large N c expansion with chiral perturbation theory (%PT), plus some additional 
assumption(s). Two most recent studies, in which such a strategy has been adopted, resulted 
in B K = 0.41(9) |Lj and B K = 0.77(5) (5) @] Q. On the other hand, the chiral quark model 
prediction of ref. [|15) favours much larger value, i.e. B^ = 1-1(2). QCD sum rules have 
also been used to compute B K , with the results dispersed over a large range of values. [] 

Clearly, these refinements of the approximate results did not help clarifying the issue of 
the correct value for B^. It should be mentioned, however, that in SM one can extract the 
range of allowed values for B%, by requiring the consistency with the S-physics constraints 
on the CKM triangle. In this way, the authors of ref. 0] found that 0.62 < Bk < 1-48 
(with 95% C.L.), thus excluding the low SM values such as B K < 0.6 (see |19| and references 
therein). 

The most direct and systematically most improvable way to compute the matrix ele- 
ment (H) in QCD is by using the lattice simulations. Over the past decade, a substantial 
progress in controlling various systematic uncertainties of this approach has been made (see 



refs. [20 , |2T| , |22| for reviews which contain a complete list of references). 

Three currently used formulations of the lattice fermion action have been employed 
in the computation of B K (or, as it is customary B^ s {2 GeV)) , and in the quenched 
approximation, they lead to the following results: 

• With the staggered fermions one can compute the matrix element (fj) by using very 
light quark masses. The renormalization of the operator Q(fi) has been performed 
perturbatively. By extrapolating to the continuum limit, JLQCD [E3j obtained 



B^ s (2 GeV) = 0.628(42) , (5) 
and thus improved (confirmed) the previous result of ref. [33], B^ s (2 GeV) = 0.62(2)(2). 

With Wilson fermions, as compared to the staggered, one cannot work with the very 
light quark masses, but the renormalization of the operator Q (which is more delicate 
due to the mixing with the operators of the same dimension) has been done non- 
perturbatively. Presently, the most refined result is the one by JLQCD ||25||. 



B™ s (2 GeV) = 0.69(7) , (6) 
fully consistent with the previous B^ (2 GeV) = 0.66(11), as obtained in ref. PB . 



1 Note that VSA and large N c cannot distinguish between Bk(h) and Bk- 

2 The result of ref. Jll| is obtained in the chiral limit in which the authors of ref. jl4| quote B\ 
0.32(6)(12). 

3 Three different QCD sum rule techniques lead to: B K = 0.5(1)(2) jl6), 0.84(8) @, 0.39(10) ||. 
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Very recently, also the domain wall fermions were employed in this computation. The 
results of RBC and CP-PACS are: 

E]f (2 GeV) = 0.538(8) [27] , 

5]f (2 GeV) = 0.575(6) . (7) 

Notice that RBC also implemented the non-perturbative renormalization in their com- 
putation m, gg. 



In summary, from the use of various lattice techniques, one concludes that the value of 
£jf(2 GeV) lies in the interval 0.53 < Bjf (2 GeV) < 0.77, or in the RGI form, 0.72 < 
B K < 1.05. Q 

In this letter, we contribute the ongoing discussion by presenting the results of our com- 
putation of the matrix element (0), for which we use the 0(a) non-perturbatively improved 
Wilson action |j3~I] , and renormalize the operator Q(n) non-perturbatively. The operator 
Q(fi) has not been improved, though. f\ In order to examine the systematic uncertainty 
arising from the use of quenched approximation, we attempted an unquenched computa- 
tion with two degenerate dynamical quarks (rif = 2) and for two choices of the dynamical 
quark masses. From the direct comparison to the quenched result, and bearing in mind the 
modest statistical quality of our unquenched data, we do not see any difference between the 
quenched and unquenched results. 

Our final numbers are: 

Bf{2 GeV) = 0.73(7)±g : g; B K = 1.01(9)1°;°! , (8) 
thus favouring the larger values for this quantity. 



2 Outline of the Lattice Computations 

In order to determine the value of the parameter B^in) in the RI-MOM scheme, we studied 
the relevant matrix element on four lattices for which the main features are summarized in 
tab. [I]. For technical details related to the production of the unquenched configurations, 
consult refs. |32], The values of the Wilson hopping parameters are listed in tab. §, to- 
gether with the corresponding masses of pseudoscalar mesons and ratios of the pseudoscalar 
to vector meson masses. In our notation the meson masses in lattice units are written in 
capital letters, whereas the small case letters are used to denote the same masses in physical 
units (e.g. M P = amp). As usual, the meson masses are extracted from the exponential 
fall-off of the two-point correlation functions. All the source operators used in this work are 

4 This range of values agrees with the one presented by L. Lcllouch at "Lattice 2000", namely 0.66 < 
B K < 1.06 |2|. 

5 As compared to ref. f26f , our action is non-perturbatively improved whereas their was tree-level im- 
proved. In addition, we do not rotate the quark fields. This rotation is unnecessary (the 4-fermion operator 
is not improved) and it introduces a statistical noise. This is why our statistical errors are smaller than 
theirs. 
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Label 


(ft 


L 3 x T) 


csw 


Statistics 


ar l {m K *)[GeV] 


Comment 


(a) 


(6.2, 


24 3 x 48) 


1.614 


200 config.'s 


2.71(8) 


Quenched 


(bl) 


(5.6, 


24 3 x 40) 





60 config.'s 


2.48(11) 


P.U. (n sea = 0.1575) 


(b2) 


(5.6, 


24 3 x 40) 





60 config.'s 


2.60(12) 


P.U. (n sea = 0.1580) 


(c) 


(6.1, 


24 3 x 40) 





60 config.'s 


2.51(7) 


Quenched 



Table 1: Characteristics of the lattices used in this work. P.U. stands for the partial unquenching 
where the value of the corresponding dynamical (sea) quark mass parameter is indicated in the 



parentheses, a 1 has been estimated from the mass of the K* meson in a way explained in ref. [ 34 j 



local. Statistical errors are estimated by using the jackknife procedure. For completeness, 
we give for each lattice the value of the critical hopping parameter as obtained by linearly 
extrapolating the square of the pseudoscalar meson mass (i.e. the valence quark mass) to 
zero. Using the same labels as in tab.0, we give: 



K 



crit 



(a) (bl) (b2) (c) 

{0.13580(2), 0.15931(6), 0.15885(7), 0.15494(5)} 



(9) 



The physical volume of the lattice (c) has been chosen to be nearly equal to the one 
used in simulations with partial unquenching (L 3 w (1.9 fm) 3 ). Moreover, the hopping 
parameters for run (c) are chosen in such a way that the masses of pseudoscalar mesons are 
very close to the ones obtained from the unquenched runs. In this way any sizable difference 
between the physical quantities obtained from the quenched and unquenched simulations 
would signal the quenching error. 

Let us now turn to the definition of the operator Q(jj) that appears in eq. (0). Since the 
Wilson term in the fermion action explicitly breaks the chiral symmetry, an extra mixing of 
the operator Q with other dimension-6 operators occurs. Consequently, in order to match 
the bare operator computed on the lattice to the continuum one, one has to subtract the 
effects of this mixing. In addition, one also needs the overall multiplicative renormalization 
constant. In a compact form, this procedure can be written as p5l |36l 



Q(fi) = Z(fi,g 2 ) \Q x + Y.^i{9l)Qi 



(10) 



where Qi on the r.h.s. denotes a bare lattice operator from the basis of 4-fermion operators, 
which we choose to be: 

Q = Qi = s^f^d S7 M d + sj^ 5 d sj^ 5 d 
Q 2 = sj^d s'j^d - s'j^d s'j^d 

Qs = sd sd — sj^d sjsd (11) 
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' " ± 


Ko 


o 


K4 


(a) 


0.1352 


0.1349 


0.1344 




M P 


0.200(2) 


0.245(2) 


0.307(1) 




Tdplvny 


0.592(23) 


0.678(16) 


0.758(9) 




(bl) 


0.1580 


0.1575 


0.1570 


0.1560 


M P 


0.230(5) 


0.268(5) 


0.303(4) 


0.366(4) 


TTlp/mv 


0.612(30) 


0.685(22) 


0.735(17) 


0.802(11) 


(b2) 


0.1580 


0.1575 


0.1570 


0.1560 


M P 


0.190(7) 


0.234(5) 


0.273(4) 


0.342(3) 


rap 1 my 


0.633(15) 


0.700(16) 


0.748(15) 


0.816(11) 


(A 


0.15363 


0.15338 


0.15311 


0.15258 


M P 


0.228(5) 


0.248(4) 


0.269(4) 


0.307(4) 


mp/m v 


0.623(12) 


0.662(11) 


0.698(10) 


0.752(8) 



Table 2: Pseudoscalar meson masses (in lattice units), and the ratio of the pseudoscalar-to-vector 
mesons consisting of two degenerate (valence) quarks of mass corresponding to the indicated value 
of the hopping parameter k%. Labels (a,bl,b2,c) refer to the lattices with characteristics given in 
tab. § 

Q 4 = sdsd + sj^d sjsd 
Q5 = sa^ u d so^d . 

To compute the subtraction (A^q)) and the renormalization (Z(fi, g^)) constants non- 
perturbatively, we use the technique proposed and elaborated in refs. [[35], |3~6| , and work 
in the RI-MOM scheme in the Landau gauge. In tab. Bl we list the non-perturbatively 
determined values for Z(n,gl) and Ai(g^) for all four lattices. The actual computation is 
performed at two values of the renormalization scale, fia = 0.7, 1.4, and for each value of the 
quark mass. For the presentation convenience, we linearly extrapolate each constant to the 
chiral limit, K q — > K CTit - This extrapolation is very smooth. In addition, we evolve Z(fia) 
to the scale /ia = 1 by using the NLO anomalous dimension coefficient (see eq. @) and by 
taking the two-loop expression for a s (fi), normalized at a s (3.41/a) by the value obtained 
from the average plaquette, as prescribed in refs. [BJ, 0. f\ At the scale fi = 1/a we make 
a comparison with the results of one-loop perturbation theory. We use the perturbative 

6 We stress that the values of Z(/j,) which are computed at two values of the renormalization scale, are 
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Label 


Method 


Z(fi = l/a,# 2 ) 


A 2 (s 2 ) 


A 3 (2 2 ) 




A 5 (S 2 ) 


(A 


NP 


0.62(1) 


-0.058(6) 


-0.018(3) 


0.016(2) 


0.003(2) 




PT 


0.68 


-0.041 


-0.008 


0.015 


0.015 


(bi) 


NP 


0.51(2) 


-0.198(16) 


-0.030(8) 


0.060(17) 


0.019(11) 




PT 


0.45 


-0.231 


-0.042 


0.084 


0.084 


(b2) 


NP 


0.52(2) 


-0.203(12) 


-0.039(8) 


0.066(18) 


0.013(10) 




PT 


0.45 


-0.230 


-0.042 


0.083 


0.083 




NP 


0.55(2) 


-0.152(8) 


-0.028(6) 


0.047(10) 


0.007(5) 




PT 


0.52 


-0.173 


-0.032 


0.063 


0.063 



Table 3: Renormalization and subtraction constants (eq. ([171)), evaluated non-perturbatively 
in the (Landau) RI-MOM scheme, extrapolated (linearly in m q ) to the chiral limit and evolved 
to fj, = 1/a by using the NLO evolution. Non-perturbative results are directly confronted to the 
perturbative values. The labels coincide with those used in tab. 0. 



expressions of ref. |39 , and combine them with the results of ref. [f|0j. In the basis (|TT| ) and 
in the RI-MOM scheme they read: 



Z(l/a, gl) = 1 - ||(3.674 - 0.742c w - 0.388c| w ) , 
Z(l/a, gl) ■ A^ 2 ) = cAo.767 - 0.795c 5H / + 0.272^) 

47T 



(12) 



where c 2 = — 11/12, c% = — 1/6, 04 = 05 = 1/3. In the numerical computations, we use the 
boosted coupling, g\ — > g\ = 6/(/3(P)) [], where the average plaquette values are computed 
in our simulations: 



(a) (bl) (b2) (c) 

(P) = {0.6136, 0.5725, 0.5735, 0.6050} . (13) 

The labels above figures are the same as in tab. [1]. From tab. |3| we see that the non- 
perturbatively computed subtraction constants Aj are not very far from their perturba- 
tively estimated values (except for the case of As((?o) whose perturbative value is generally 

consistent with the NLO anomalous dimension (J3J) - The difference that arise from the evolution from the 
two scales to /1 = 1/a is smaller than the statistical error and it has been included in the errors quoted in 
tab. I 

7 In the improved case, we use csw = 1, to be consistent with the 1-loop order perturbation theory. 
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much larger than the non-perturbative one). The value of the non-perturbatively computed 
renormalization constant differs from the perturbative one by (10 -j- 15)%. 



2.1 Computation of the Bk{h) parameter 



The strategy to compute the bag parameter is well known [gT| and we only briefly recall it. 
The asymptotic behavior of the two- and three-point correlation function is given by the 
following expressions: 



C#(a;t) = E e-^V(f,t)jt(o,o)> 



z 



2E, 



J_ e -Ejt 



xy 



-t Pll t P2 »0 VZ P _ e _ Ep{qltpi { p mQ ^ )lp{m ^l^ p -E P (v)t P2 



2E P (q) 



2E P (p) 



(14) 



where the operator Q(fjL) is given by eq. (PH), \fZp = (0|P 5 |P) (P 5 = iq^q), p, q are 
momenta given to the interacting pseudoscalar mesons. In all our simulations, we fix tp 2 = 
12, while tp l is free. The elimination of the exponentials in the three-point function is 
achieved by considering the ratio 



R(p, q- //) 



C (3) (p,g;t Pl ,tp 2 ;/i) 



Z\ Cp P {p;tp 2 ) Cpp(q;tp 1 



(2) 



(15) 



We also divided by the axial current renormalization constant which we computed non- 
perturbatively by using the method of ref. [32]: R 



(a) (bl) (b2) (c) 

Z A (gl) = {0.80(1), 0.81(1), 0.80(1), 0.82(1)} . (16) 



After a careful inspection of the ratio (|T5| ) for several different values of p and q on the 
lattice (a), we select the following non-equivalent combinations: 

p= (0,0,0) & q = {(0,0,0), (1,0,0)}, 



p= (1,0,0) & q = {(1,0,0), (0,1,0)}, (17) 

where each component is given in units of (2ir/L). For the lattices (bl), (b2) and (c), 
we consider only the first two combinations. This is because the temporal axis of these 
lattices is rather small (T = 40) to accommodate a clear signal when the larger momenta 
are given to the 'kaons'. As usual, whenever possible, we average over equivalent kinematical 
configurations. 

8 Recent works in which Za and other bilinear quark operators' renormalization constants are computed 
by using the method of ref. H), can be found in rcfs. Q fl|, In refs. |§ 0, has also been 
computed by using the hadronic Ward identities. The results of the two methods are consistent. 
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Since the Wilson term in the fermion action breaks chiral symmetry, in the chiral ex- 
pansion of the matrix element (0), one also allow for a free term and the one proportional 
to the quark mass (i.e. to m 2 K ): 



(K (q)\Q(n)\K°(p)) = a + (3 m\ + 7 (p • q) + • • • 



where the dots stand for higher order terms in expansion. The coefficients a and (3 in 
eq. (0) are expected to be consistent with zero if the chiral behaviour of the operator Q(fi) 
is correct. 

To investigate this point we follow ref. [[IT] and rewrite eq. (|TJ) to the first order, in the 
following form: 



(19) 



where the suitable choice for the fit with lattice data is 

> Un(.r)^-.(U)> 2 



X 



£<p 5 (*)p 5 (o)> 



and 



Y 



p-q 



X 



m 



(20) 



K 



It is now clear that the division by Z\ in the ratio (|15|) provides the renormalization of 
the axial current (A = 57075(7) in the quantity X, which allows then to identify the term 
proportional to Y as f\{j> • q), where fx is the (renormalized) kaon decay constant defined 
from (0|^|#°(p)) = if KPft . 

Before discussing the physical results, we give in tab. |] and |5] the explicit values of the 
quantity Y and the ratio R(p,q;fi). We stress that the renormalization and subtraction 
constants used in this computation are those obtained for each value of n q separately. We 
have also computed R(p, q; fi) by using the renormalization and subtraction constants ex- 
trapolated to the chiral limit. Since the extrapolation of these constants to the chiral limit is 
very smooth, 9 the net effect on our results is negligibly small (less than 1%). In computing 
(p ■ q) we used the latticized energy-momentum relation 



sinh 



sinh 



Mr 



(21) 



which describes well our data 34 . 



3 Physical results 

Our central result is obtained from the fit of our lattice data (a), listed in tab. |], to the 
form given in eq. (|T9"D. This fit is illustrated in fig. [I] and the results are: 



/ia = 0.7 



a 



-0.017(10); (3 = 0.10(10); Bf~ MOM {fi) = 7 = 0.713(67), 



3 The corresponding plots can be obtained from the authors. 
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(a) 



P, Q 


(0,0,0), (0,0,0) 


(0,0,0), (1,0,0) 




y = x 


R(p,q) 


y 


R(p,q) 


Ki 
^2 
% 


0.068(7) 
0.102(8) 
0.166(10) 


0.040(6) 
0.068(8) 
0.119(12) 


0.111(11) 
0.149(12) 
0.217(13) 


0.068(10) 
0.102(10) 
0.159(12) 




p, q 


(1,0,0), (0,1,0) 


(1,0,0), (1,0,0) 




y 


R(p,q) 


y 


R(p,q) 


«1 
^2 
«3 


0.182(17) 
0.216(17) 
0.283(17) 


0.111(17) 
0.150(18) 
0.205(18) 


0.066(7) 
0.100(8) 
0.163(10) 


0.033(17) 
0.063(16) 
0.108(14) 



Table 4: Numerical values of the quantity Y and of the ratio R(p, q; (J,), computed in the RI-MOM 
scheme at fia = 0.7 (fi = 1.9 GeV) for the lattice (a). The results are presented for each value of 
the hopping parameter Ki specified in tab. [|, and for four combinations of momenta p and q. 



H<i = lA : a = -0.009(9); (3 = 0.10(10); Bg- MOM (n) = 7 = 0.706(65). (22) 

To convert J B| I_MOM (/i) to the RGI form and also to the MS (NDR) scheme, in which the 
lattice results are usually presented, we use eq. flU). At this point we should decide which 
coupling to use. One choice is to take the value ai n/_ °' ) (l/a) = 0.178 extracted from the 
average plaquette (P), by using the recipe of ref. ||38|| . Another (more phenomenological) 
choice is to relate the quenched S^ I ~ MOM (/i) to the physical (unquenched) world by using the 
physical ai" /_4 ' ) (l/a) = 0.266, which is obtained by the NLO evolution of a s {M z ) = 0.118. 
The difference between the two conversions is less than 1%. In this letter we use ai n/_ °^(l/a), 
which leads to 

5jf (2 GeV) = 0.730(68); B K = 1.009(94) . (23) 

As mentioned before, the results obtained at two different scales (|22|) are fully consistent with 
the first two coefficient of the anomalous dimension obtained in perturbation theory. The 
numbers given in eq. (|23|) are the average of the conversions i?^ I_MOM (/i) — > B^ s (2 GeV) 
from two different values of the scale /1. This is the main result of this letter to which we 
should add the systematic uncertainties which we now discuss. 

- Our pseudoscalar mesons are composed of degenerated quarks. The quenched 
suggests that the value of Bk{^) may increase by (3 4- 5)% if we work with non- 
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(bl) 



p, q 


(0,0,0), (0,0,0) 


(0,0,0), (1,0,0) 


Ki 


Y = X 


R(p,q) 


Y 


R(p,q) 


Ki 


0.078(3) 


0.070(26) 


0.117(6) 


0.101(34) 




0.108(3) 


0.091(27) 


0.151(5) 


0.120(34) 


K 3 


0.113(3) 


0.113(28) 


0.183(5) 


0.143(30) 


K4 


0.200(4) 


0.162(27) 


0.244(5) 


0.189(33) 


(b2) 


p, q 


(0,0,0), (0,0,0) 


(0,0,0), (1,0,0) 


Ki 


Y = X 


R(p,q) 


Y 


R(p,q) 


K X 


0.060(5) 


-0.002(22) 


0.101(10) 


0.013(40) 


K 2 


0.090(5) 


0.024(22) 


0.134(8) 


0.051(35) 


K 3 


0.121(5) 


0.051(24) 


0.166(7) 


0.085(34) 


K4 


0.183(4) 


0.108(28) 


0.229(5) 


0.148(34) 


(c) 


p, q 


(0,0,0), (0,0,0) 


(0,0,0), (1,0,0) 


K{ 


Y = X 


R(p,q) 


Y 


R(p,q) 


Ki 


0.083(14) 


0.056(28) 


0.125(16) 


0.083(27) 


«2 


0.100(15) 


0.069(26) 


0.143(16) 


0.097(26) 


«3 


0.118(16) 


0.083(24) 


0.163(17) 


0.114(25) 


Ki 


0.155(17) 


0.113(21) 


0.202(20) 


0.148(24) 



Table 5: The same as in tab. [| but for lattices (bl), (b2) and (c) in which we consider two 
combinations of the momenta p and q. 
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Figure 1: Fit of our lattice data (a) to eq. (|T^), The shown values of R(p,q\fi) are those corre- 
sponding to fia = 0.7, computed in the (Landau) RI-MOM scheme. The fitting curves (dashed 
lines) correspond to three different values of X(ni) in eq. (|T^), 



degenerate quarks f48fl . To be on a safe side, we include +5% to the systematic 
uncertainty. 

- We tried to modify the expansion fll9D , by including the term X\og(mp/ fp) which is 
the one that shows up in both the standard and the quenched xPT . We also tried 
to modify (0) by adding a term proportional to X 2 . Any of these two modifications 
does not change the value of Bk{^)- 

- If instead of the procedure explained above, we use the renormalization and subtraction 
constants extrapolated to the chiral limit before doing the fits to eq. fll9|), the final 
value for Bk(h) gets lower by less than a percent. To be conservative, we take ±1% 
which would also account for the tiny difference arising from the choice of a s (/j) in 
the conversion formulae. 

- Since we work at fixed value of the lattice spacing it is difficult to estimate the dis- 
cretization effects. Even though the action we use is free of 0(a) effects, our operator 
Q(f-i) is unimproved. To get some idea about the remaining 0(a) uncertainty, we re- 
peated the analysis by only improving the axial current, J^ mpr ^ = Aq + CAd^P^ (with 
the value of ca = —0.04 pi|). This further modifies the final Bk(h) by +1%. R 



10 It can also be deduced from ref. [^5| that the difference between Bk as obtained at /3 > 6.1 and the one 
extrapolated to the continuum limit is very small. 
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- Finally, one wonders about the size of the systematic uncertainty that comes from 
the use of quenched approximation. A comparison between the chiral expansion for 
the Bk parameter in the quenched and in the standard xPT, indicates that the effect 
of quenching might be negligible. In the next subsection, we will show that our 
unquenched data indeed tend to confirm that expectation. 

By simply summing up all thee above (small) uncertainties, we end up with: 

Bf{2 GeV) = 0.73(7)±8:gg; B K = 1.01(9)1°$ , (24) 
which is our final result. 

3 . 1 Unquenching 

We now discuss the results of the analysis of our unquenched data. To our knowledge, this 
is the first attempt to compute the unquenched Bk((i), by using Wilson fermions. In fact, 
our results are only partially unquenched (P.U.) which means that a valence quark mass is 
different from the dynamical one. Our both sets of data are obtained with nf — 2 ((bl) 
and (b2) in tab. p]). As compared to the physical strange quark mass, the dynamical (sea) 
quark masses are, mf^ /m s ~ 0.75, and mf^ jm s ~ 0.50. The results of the fit to eq. (|T^) 
are given in tab. || and illustrated in fig. (H). As in the previous subsection, we convert the 



Label 


fia 


a 


P 


7 = Bf- MOM (fi) 


(bl) 


0.7 
1.4 


0.01(3) 
-0.01(3) 


0.06(19) 
0.02(16) 


0.67(20) 
0.65(17) 


(b2) 


0.7 
1.4 


-0.06(3) 
-0.08(3) 


0.28(24) 
0.21(24) 


0.66(26) 
0.65(26) 


M 


0.7 
1.4 


-0.01(3) 
-0.02(3) 


0.13(19) 
0.09(19) 


0.69(18) 
0.69(18) 



Table 6: Result of the fit of our data obtained from the computation on the lattices (bl), (b2) 
and (c), to the form given in eq. |lS|) . 

obtained results to the RGI form and to the MS scheme. We use ai n/_2 ^(l/a) = 0.223 and 
0.222 for lattices (bl) and (b2) respectively, and a^ f (1/a) = 0.187 for the lattice (c), 
which as in the previous subsection are obtained from the average plaquettes. We finally 
have: 

(bl) B^{2 GeV) = 0.67(19); B K = 0.94(26) , 
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Figure 2: The same as for fig. [T], but for the partially unquenched data (bl). 



(b2) B^ s (2 GeV) = 0.67(26); B K = 0.94(37) , 
(c) B K m {2 GeV) = 0.70(18); B K = 0.97(25) . (25) 

The statistical errors in these results are large. Systematics due to the lattice geometry is 
the same for both sets of data (quenched and unquenched), so that the difference would 
indicate the quenching errors. 

First of all, from the above numbers ((bl) and (b2)) we do not observe any difference 
in Bjc as the value of the sea quark mass changes. In principle, to get a fully unquenched 
result, one should work with more than two values of the dynamical quark mass, and then 
extrapolate Bk in such a way that the dynamical quark mass corresponds to the physical 
up/down quark. Since we do not see any change in B K (/u,) as the value of the sea quark 
mass change, we can not distinguish partially from the fully unquenched value for B K (/i). 

Even more, we do not see any significant effect due to the unquenching, i.e. by switching 
off the quark loops from the simulations (bl) and (b2), we get the results from the run (c), 
which totally agrees with the unquenched ones. This observation is on the same line with 
the conclusions drawn from the comparison of the standard and the (partially) quenched 
%PT [[|^, [50 1 • In re f- |51 1) however, by using the staggered fermions, an enhancement of Bk 



of the order of (5 ± 2)% due to the unquenching has been reported. Obviously, with a poor 
statistical quality of our data, we are unable to confirm this result. 

11 In fact, the enhancement in Bk due to the unquenching ~ 5%, as claimed by the authors of ref. p^ , 
comes from their simulations with nf = 4, while their data with rif — 2 almost coincide with the quenched 
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In conclusion, we do not add any systematic uncertainty to our result due to the 
use of quenched approximation. 



4 Conclusion and prospects 

In this letter, we computed the Bk parameter by using 0(a) improved Wilson action but 
without improving the corresponding four fermion operator. From the non-perturbatively 
renormalized data, we obtained B^ s (2 GeV) = 0.73(7)io]oi) which confirms the tendency 
of larger values for this quantity when working with the Wilson fermions [ G5L |39], p^ ]. 
In order to examine the effect of quenching, we computed Bk with dynamical fermions 
(rif = 2). Even though the statistical quality of our unquenched results is poor, we do 
not see any significant deviation of the parameter Bk as compared to the one obtained in 
the quenched simulation. To quantify the small quenching errors (if any!) it is obviously 
important to increase the statistics of the unquenched data. 

The main systematic uncertainty, which we could not estimate, is expected to come from 
the computation of the bare operator, i.e. from the subtractions of the effects of mixing 
with operators of the same dimension (peculiarity of the Wilson lattice regularization). This 
can be completely avoided if one follows the proposal of ref. ||53|| , namely by using the Ward 
identity to relate Q(f-i) = VV + AA to the parity violating operator = V A + AV 

for which no subtraction is needed. Alternatively, one can use the twisted mass fermions, 



as proposed in ref. [54]. It will be very interesting to see if the results for Bk without 



subtractions will confirm the values obtained in this letter. 



Acknowledgements 

We thank T\T Collaboration for providing us with the unquenched gauge field configura- 
tions. We are grateful to V. Gimenez, V. Lubicz and G. Martinelli for many discussions, 
and to Ph. Boucaud, L. Giusti and A. Vladikas for their comments on the manuscript. D.B. 
also thanks INFN for the financial support. 



References 

[1] D. E. Groom et al, Eur. Phys. J. C15 (2000) 1. 

[2] A. J. Buras,in "Probing the Standard Model of particle interactions", Les Houches, 
Qhep-ph/980647lQ . 

[3] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125, 
l jriep-ph/9512380p . 

[4] L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945. 
ones. 



14 



[5] A. J. Buras, M. E. Lautenbacher and G. Ostermaier, Phys. Rev. D 50 (1994) 3433, 
[ |hep-ph/9403384| ] . 



[6] 
[7] 

[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 
[17 
[18 

[19 

[20 
[21 

[22 

[23 

[24 

[25 



M. K. Gaillard and B. W. Lee, Phys. Rev. Lett. 33 (1974) 108. 

M. Ciuchini, E. Franco, V. Lubicz, G. Martinelli, I. Scimemi and L. Silvestrini, Nucl. 
Phys. B523 (1998) 501, [ |hep-ph/97lT402| ] . 



A. J. Buras, M. Misiak and J. Urban, Nucl. Phys. B586 (2000) 397, Qhep-ph/0005183| ] . 
G. Altarelli, G. Curci, G. Martinelli and S. Petrarca, Nucl. Phys. B187 (1981) 461. 
A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B347 (1990) 491. 
J. F. Donoghue, E. Golowich and B. R. Holstein, Phys. Lett. B119 (1982) 412. 
W. A. Bardeen, A. J. Buras and J. M. Gerard, Phys. Lett. B211 (1988) 343. 



S. Peris and E. de Rafael, Phys. Lett. B490 (2000) 213, [ |hep-ph/0006l"4l ] . 



J. Bijnens and J. Prades, JHEP0001 (2000) 002, [ |hep-ph/99lT39l ] . 

S. Bertolini, J. O. Eeg, M. Fabbrichesi and E. I. Lashin, Nucl. Phys. B514 (1998) 63, 
[ |hep-ph/9705244r j . 



N. Bilic, C. A. Dominguez and B. Guberina, Z. Phys. C39 (1988) 351. 
L. J. Reinders and S. Yazaki, Nucl. Phys. B288 (1987) 789. 

J. Prades, C. A. Dominguez, J. A. Penarrocha, A. Pich and E. de Rafael, Z. Phys. C51 
(1991) 287. 

M. Ciuchini, G. D'Agostini, E. Franco, V. Lubicz, G. Martinelli, F. Parodi, 
Uhep-ph/0012308p . 



S. R. Sharpe, Nucl. Phys. Proc. Suppl. 53 (1997) 181, [ |hep-lat/9609029p 



L. Lellouch, talk given at "Lattice 2000", Bangalore, 17-22 August, 2000, 
L |riep-lat/0011088| ]. 



V. Lubicz, talk given at "20 th Physics in Collision Conference", Lisbon, Portugal, 29 
Jun - 1 Jul 2000, Qhep-ph/00f0l7lt ] . 

S. Aoki et al. [JLQCD Collaboration], Phys. Rev. Lett. 80 (1998) 5271, 
[ |hep-lat/9710073| ] . 

G. Kilcup, R. Gupta and S. R. Sharpe, Phys. Rev. D 57 (1998) 1654' 
L [hep-lat/9707006| ] . 

S. Aoki et al. [JLQCD Collaboration], Phys. Rev. D 60 (1999) 034511, 
Qhep-lat/9901018| ]. 

15 



[26] L. Conti, A. Donini, V. Gimenez, G. Martinelli, M. Talevi and A. Vladikas, Phys. Lett. 
B421 (1998) 273, Qhep-lat/9711053p . 



[27; 

[28 
[29 
[30 
[31 

[32 

[33 

[34 
[35 
[36 

[37 
[38 

[39 

[40 
[41 
[42 

[43 
[44; 



T. Blum et al. [RBC collaboration], L [hep-lat/0011042[ ] . 

A. Ali Khan et al. [CP-PACS Collaboration], [ |hep-lat/0010079p . 

C. Dawson [RBC Collaboration], Qhep-lat/0011036| ] . 
Y. Zhestkov, [ }iep-lat/0011002| ] . 

M. Luscher, S. Sint, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B491 (1997) 323, 
l ]hep-lat/9609035| ] . 

T. Lippert et al. [SESAM Collaboration], Nucl. Phys. Proc. Suppl. 63 (1998) 946, 
Qhep-lat/9712020| ] . 

N. Eicker et al. [SESAM collaboration], Phys. Rev. D 59 (1999) 014509, 
Qhep-lat/9806027| ] . 

D. Becirevic et al, [ |hep-lat/9809129P . 

M. Crisafulli et al, Phys. Lett. B369 (1996) 325, [ |hep-lat/9509029| ] . 

A. Donini, V. Gimenez, G. Martinelli, M. Talevi and A. Vladikas, Eur. Phys. J. CIO 
(1999) 121, [ |hep-lat/9902030| ] . 



G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 48 (1993) 2250, [ ]hep-lat/9209022D . 

C. T. Davies, K. Hornbostel, G. P. Lepage, P. McCallum, J. Shigemitsu and J. H. Sloan, 
Phys. Rev. D 56 (1997) 2755, Qiep-Iat/970301CD ■ 

R. Gupta, T. Bhattacharya and S. Sharpe, Phys. Rev. D 55 (1997) 4036, 
Qhep-lat/9611023| ]. 

S. Capitani, et al, Qhep-lat/971 10071 ] . 

M. B. Gavela, et al, Nucl. Phys. B306 (1988) 677. 

G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa and A. Vladikas, Nucl. Phys. B445 
(1995) 81, Uhep-lat/9411010p . 



M. Gockeler et al, Nucl. Phys. B544 (1999) 699, Qhep-lat/9807044f ] . 

V. Gimenez, L. Giusti, F. Rapuano and M. Talevi, Nucl. Phys. B531 (1998) 429, 
Qhep-lat/9806006| ] . 



[45] D. Becirevic, V. Lubicz, G. Martinelli and M. Testa, Nucl. Phys. Proc. Suppl. 83 (2000) 
863, [[hep-lat/9909039|] . 



16 



[46] M. Luscher, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B491 (1997) 344, 
Qhep-lat/9611015P . 

[47] T. Bhattacharya, R. Gupta, W. Lee and S. Sharpe, [ Jhep-lat/0009038| ] . 

[48] S. R. Sharpe, in "CP violation and the limits of the Standard Model", TASI 94, (World 
Scientific), [ [hep-pri/9412243; ] . 

[49] S. R. Sharpe, Phys. Rev. D 46 (1992) 3146, [ |hep-lat/920502C| ] . 

[50] M. F. Golterman and K. Leung, Phys. Rev. D 57 (1998) 5703, Qhep-lat/9711033p . 

[51] G. Kilcup, D. Pekurovsky and L. Venkataraman, Nucl. Phys. Proc. Suppl. 53 (1997) 
345, Qhep-lat/9609006p . 

[52] L. Lellouch and C. J. Lin [(UKQCD Collaboration)], Nucl. Phys. Proc. Suppl. 73 (1999) 
312, [ piep-lat/9809142[ J . 

[53] D. Becirevic, P. Boucaud, V. Gimenez, V. Lubicz, G. Martinelli, J. Micheli and M. Pa- 
pinutto, Phys. Lett. B487 (2000) 74, [ ^iep-lat/0005013| ] . 

[54] R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz, Nucl. Phys. Proc. Suppl. 83 (2000) 
941, [[hep-lat/9909003|] . 



17 



